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Abstract 

A number of experimental measurements of the Casimir force have observed a logarithmic dis- 
tance variation of the voltage that minimizes electrostatic force between the plates in a sphere-plane 
geometry. We show that this variation can be simply understood from a geometric averaging of 
surface potential patches together with the Proximity Force Approximation. 
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A number of experimental measurements of the Casimir force have observed a dis- 
nce variat: 
potential, jl- 



tance variation of the voltage applied between the plates that minimizes the electrostatic 
-|4| This distance variation is of the approximate form 



V m (d) = a + b\ogd (1) 

over a range of distances d spanning up to nearly two orders of magnitude. We have shown 
numerically that a variation in the minimizing voltage can result from the geometrical aver- 
aging of patch potentials on the plate surfaces, and have developed a heuristic explanation 
of the effect, as described in js|. However the logci form of the variation was not explained 
or derived in j^]. Here we provide added details showing that in the case of a spherical 
surface/plane surface geometry, the logo? variation arises quite naturally. 

The electrostatic force between the plates is minimized when the free energy is minimized. 
For a plate with a spherical surface curvature R together with a flat plate, both of diameter 
2R m , the Proximity Force Approximation (PFA) works well for distances d < R? m j2R\ in 



addition, the characteristic radius of the patches r must satisfy r > v2Rd, in which case 
we can consider each patch as only interacting with its own image in the opposite plate. 
If this latter criterion is not met, the effects calculated here will be reduced in magnitude, 
however the general conclusions are otherwise unaltered. 

For random patches on the surfaces, the electrostatic free energy is 

where V a is a voltage applied between the plates, and V p (r, <fi) describes the random voltage 
patches on the plates' surfaces. The value of V a that minimizes the force (or electrostatic free 
energy) is referred to as the minimizing potential, and is often called the contact potential. 
V m is found by taking the derivative with respect to V a : 

dU A n f Rm f 2n nr j.^( r / R ) dr 

W a 

Because the only <p dependence is in V p (r, (ft), the angular integral of this term can be replaced 
with V p (r), with the understanding that this replacement represents a sum of the individual 
patches, of typical radius ro, intersected by a circle of radius r. For homogeneous random 
patches, this number is roughly (2nr) / (2r ) , for r > r , or 1 for r < r . Then if RMS 
magnitude of the patches is V with zero average, this sum will have magnitude proportional 



= e oR £ m J\v a - VArA)) ^^ for V a = V m . (3) 



to the square root of the number of patches in the sum times ±Vo, or 



2tt 



{V m - V p {r, (j)))d4> = 2vr [V m + V p {r)\ (4) 
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A0 [9{t - r ) - 9{r - R m )} \l^r + ^ W) - 0(r - r )] 



V (5) 



where 6(r — tq) is the Heaviside step function, A0 = 27rro/r results from dip in the integral, 
and a single patch near the center of the plates is included. (The relative sign of these 



two terms on the r.h.s. is random.) Thus, V p (r) oc 1/yr/ro for r ^> tq. Although we 
will not directly use this result in the subsequent discussion, it is important to note that 
< l^p(0)| in the case of homogenous random patches. 
Equation ((3]) can be integrated by parts, yielding 

/Rm Ay ( r .\ 

\og(d + r 2 /2R)^r 1 dr= (6) 

= (V m - V p (R m )) log(d + R 2 J2R) - (V m - V p (0)) \ogd- Q(d) = (7) 
where we have introduced a new function Q(d). Then 

_ V p (R m ) log(d + R 2 J2R) - V p (0) log d + Q(d) 
m[) ~ \og(d + Rl/2R)-logd • 18 > 

There are three cases to consider. 
Case 1: Close range: d « R 2 m /2R } \ \ogd\ » \ hgR 2 n /2R\ and d « r%/2R 

In this limit, one single patch at or near the center dominates in the determination of V m , 
and dV p (r)/dr when r 2 /2R < d in the integral defining Q(d). The logd terms have the 
largest magnitude, and thus we can neglect all terms not multiplied by log d. In this case, 

V m (d) = V p {0). (9) 

Case 2: Intermediate range: d « Rf n /2R, \ log c?| < | logi?^/2i?| and d < r^/2R 
In this case, 

\og(d + R 2 m /2r)Klog(R 2 j2R) (10) 

so 

Q(d) = Q log(R 2 j2R) (11) 

because in this limit Q(d) is nearly independent of d. This can be seen in Eq. (5) where the 
log(d + r 2 /2R) factor contributes a d dependence only when r is small, but in that limit, 



dV p (r)/dr 0. Therefore, 



T/ „ V P (R m ) log(R 2 m /2R) - VM logd + Q log(R 2 m /2R) 

Vm{d) log(i£/2tf)-logd • (12) 



The denominator can be Taylor expanded and we therefore arrive at 

logrf 



V m (d) 



V(R v V P (0)logd 



1 + 



a + b\ogd (13) 



M^/2i2). 

where a and 6 do not depend on d, and where terms only first order in log dj \og(R% n /2R) 
are retained. In addition, there can be a contribution to a from external circuit contact 
potentials. 

It should be noted that a single patch at large r will generate a non-zero Qo. Depending on 
the size and magnitude of the single patch, it might dominate the log d distance dependence 
of V m . 

Case 3: d > R 2 J2R and d ~ r\j2R (so that the PFA remains valid) 

For this case, it is easiest to go back to Eq. (3) and expand the denominator. We have 

= / / {y m -V p (r J <f ) ))r{l-r 2 /2Rd)drd<t ) ^ / / (V m - V p (r, (f)))rdrd(j) (14) 
Jo Jo Jo Jo 

so that 

which is the surface average of V p (r, 0). 
Discussion 

It is easy to understand the first and third cases where d — > and d — > oo. In the first 
case, a single patch dominates the electrostatic force, and its potential determines V m . In the 
latter, V m is simply the average surface potential, as the variation is the distance between 
the surfaces due to the curvature is very small compared to d. 

The intermediate case is slightly more difficult to understand, but arises from, with 
increasing distance, the loss in dominance (in magnitude) of log d, compared to log(d + 
R 2 m /2R). 

Numerical calculations based on Eq. (3), using random patches specified on a surface, 
are straightforward and fully support the essential conclusions presented above. As an aside, 
there is a remaining question as to where it is the energy or force that must be minimized. 
Numerical calculations based on the minimization of the force produces no statistically 
significant differences compared to the energy minimization, as might be expected. 



Of course, this is a very simplistic model, particularly in the assumptions that the patch 
potentials are randomly ±Vb, that all patches have the same radius and are circular, and that 
the radius of the positive patches is equal to that of the negative patches. These assumptions 
are not essential to the derivation of Eq. (8), which is the principal result reported here. 
Further refinements will likely not significantly change the subsequent conclusions presented 
here. For example, as discussed already, a single patch at large r will generate a non-zero 
Qq. If this patch has a large area and/or large potential {e.g., a charged speck of dust or a 
scratch), it could easily be the dominant contribution to the logrf distance dependence. For 
such a patch, V m {d) will only vary slowly with a relative translational repositioning of the 
plates. 
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